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Abstract 

First principle calculations have been performed to understand the experimentally observed size 
sensitive variations in the characteristics of heat capacities of gallium clusters [G. A. Breaux et. 
al. J. Am. Chem. Soc, 126, 8628 (2004)]. It was reported that while some clusters exhibit a 
clear solid like to liquid like transition others exhibit a continuous transition with no peak in the 
heat capacity curve. In addition, the clusters also exhibit a variation of about 300 K (500-800 K) 
in the melting temperature across the size range of 20 to 46. In the present work we correlate the 
observed finite temperature properties to its geometry and nature of bonding in the ground state. 
We demonstrate that the local order (i.e., island of atoms bonded with similar strength) in the 
ground state geometry is responsible for the variation in the shape of the heat capacity curve. We 
attribute the higher melting temperature of clusters to the presence of distinct core and strong 
covalent bonds between the core and surface atoms. 

PACS numbers: 61.46.Bc, 36.40.Mr, 36.40.-c, 36.40.Cg 



I. INTRODUCTION 



During the past decade or so, a number of systematic calorimetric measurements probing 
the finite temperature behavior of unsupported free clusters have yielded unexpected and 
rich physics.— >2»M One of the first reports is on the measurement of melting temperatures 
(T m ) of simple metal clusters of sodium in the size range of 55 to 357 atoms.— The measured 
melting temperatures were found to be an irregular function of the cluster sizes with melting 
peaks showing no clear correlations with the observed magic numbers in sodium. 2 The 
second measurement was reported by Jarrold and co-workers which brought out additional 
important aspects in the melting of finite size clusters.—"^ The systems investigated are free 
clusters of gallium in the size range of 30 to 55.— The measured heat capacities of the clusters 
in the above range revealed three very interesting features: (1) Higher than bulk melting 
temperatures (T m [buik] — 303 K) in all the clusters (also seen in two smaller clusters viz., 
Gai7 and Ga2o) in direct contradiction with the accepted paradigm, viz., the reduction 
in the melting temperatures with the size.— (2) The size sensitive behavior of the shape 
of the heat capacities where addition of even one atom results in a dramatic change of 
shape, prompting some of the clusters to be called as "Magic Melters".— This means that 
while some clusters do undergo a conventional and clear melting transition, others undergo 
a near continuous transition making it very difficult to identify any meaningful transition 
temperature. (3) A variation of about 300 K (500-800 K) in the T m across the sizes between 
20 to 46. For example, Ga2o, Ga46, and Ga47 exhibit a highest melting temperature of about 
800 K while the intermediate clusters such as Ga3i, Ga33, Ga37, Ga4i, etc., exhibit a lower 
melting temperature of the order of 500-600 K. Thus, the size sensitive variations in melting 
temperature are not monotonic in nature. Such variations in the melting temperature of 
clusters are also observed experimentally in the case of aluminum clusters.— 

An explanation and understanding of some of the unexpected and puzzling experimental 
observations noted above warranted detailed ab initio Molecular Dynamical (MD) simula- 
tions based on Density Functional Theory (DFT). A clear signal for needing such an approach 
came from detailed investigations using classical inter atomic potentials, which failed to re- 
produce the data on sodium, even qualitatively— On the other hand, even early DFT based 
investigations, limited to short simulation times, were successful in providing some insights 
into the phenomena of melting.— Recently, works of Chacko et. al.,— and Aguado and Lopea^ 
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over the entire size range have reported excellent agreement with the experimental melting 
temperatures of these clusters. These calculations indicate that the geometry plays a more 
direct and significant role with the electronic structure influencing it in a more subtle way. 
Indeed it must be emphasized that the nature of bonding and the energetics turn out to 
be two crucial ingredients, making an explicit quantum mechanical treatment of electrons 
mandatory. 

In two short communications, we have addressed the issues of higher than bulk melting 
point — and size sensitive nature observed in gallium clusters.— The predominantly covalent 
nature of bonding in the gallium clusters in contrast with predominant metallic bonding in 
bulk-Ga has been shown to be responsible for the higher than bulk melting point of clusters. 
These reports are limited to Gai3 and Ga.n which was later on extended to Ga3o and Ga3i. 
By examining the geometry and the finite temperature behavior of Ga3o and Ga3i, the pair 
displaying maximum size sensitivity, we explained that the dramatic change in their heat 
capacities has geometric origin. 11 Specifically we established a direct correlation between 
the nature of the ground state and the observed heat capacity. An "ordered" cluster is 
expected to display a well characterized melting transition showing an identifiable, albeit 
broad peak. On the other hand a completely "disordered" cluster will undergo a continuous 
transition with a very broad heat capacity. Quite clearly such a description is qualitative 
and needs further quantifications. It also raises questions about the universality of such a 
phenomenon. The earlier work is based on two specific clusters, viz., Ga3o and Ga3i (as well 
as Gai 7 and Ga2o)^ where such an effect was discussed and did not address the issue of 
observed systematic behavior over entire range and the shift in T m to higher value in case 
of Ga 2 o and Ga 46 . 

The present work is an attempt to understand and evolve a complete picture of the 
experimental observations over the reported size range. Evidently, complete thermodynamic 
calculations based on an ab initio MD are prohibitively expensive. However, our earlier work 
clearly brings out the two main ingredients which are essential and perhaps sufficient for 
such an understanding, namely the analysis of ground state geometry and the nature of 
bonding in the clusters. In the present work, we have carried out such an analysis for 
the gallium clusters in the size range of 20 to 55. It may be mentioned that many of 
the geometries do not display any obvious rotational symmetry. This is especially true 
of large clusters. Therefore, it is not straightforward to characterize the nature of order 
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and disorder using some quantitative measure. In this sense the situation is some what 
unsatisfactory. However, considerable insight can be gained by analyzing the shape of the 
cluster, distribution of shortest bond lengths (in other words strongest covalent bonds) and 
pair correlation function. Another extremely useful quantity is the Electron Localization 
Function (ELF) which provides a semi quantitative picture of the nature of bonding and the 
connectivity of atoms. In view of the fact that we have made use of ELF extensively, this 
function will be discussed in detail later. 

It is most profitable to put the present work on gallium clusters in the context of known 
structural trends in small gallium clusters and bulk properties. The evolutionary trends in 
the geometries of small clusters ( N < 26) have been investigated by Song and Cao.— In 
smaller clusters it is relatively easy to discern the symmetry. For example, the ground state 
geometry of Gai3 is decahedron. The next two structures in size are obtained by capping 
one of the square faces by one and two atoms, respectively. This distorts the structure 
and disturbs the symmetry. By the time four atoms are added, viz., Gai7, the structure 
is completely distorted. This evolutionary trend changes after Gai7 and a structure with 
identifiable ordered planes develops e.g., as seen in Ga2o- Such an evolutionary trend in 
the geometries from well ordered one like Ga 13 to another partially ordered one like Ga2o 
through a disordered one like Gai7 is responsible for the size sensitive nature. Since there 
is nothing specific about gallium per se, in so far as the evolutionary trends are concerned, 
the size sensitive nature of the heat capacity should be an universal phenomenon. Indeed, 
our calculation on clusters of Na n (n= 55 and 58),— Au n (n=19 and 20), 15 and smaller 
clusters of Ga„ (n=17 and 20) 12 supports this conjecture. In light of the above discussion, 
the availability of experimental heat capacity data on gallium for all the sizes within a given 
range presents an interesting opportunity to carry out detailed investigations of their ground 
state properties. 

In the bulk phase, gallium has been described as a metallic molecular solid. Gong et. 
al.,J£ have calculated the band structure of a-Ga and have shown that molecular character 
and metallic conduction coexist. a-Ga is a face centered orthorhombic solid with 8 atoms in 
the cubic unit cell. It has only one nearest neighbor at a distance of 2.44 A. The next three 
shells contain two atoms each at distances of 2.71 A, 2.74 A, and 2.83 A, respectively. The 
structure can be regarded as strongly buckled planes connected by the shortest bonds. Fig. [T] 
shows the atomic arrangement from this perspective. The short bonds across the planes 
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FIG. 1: Unit cell of a-Ga. It shows two buckled planes. The dark line joining the atoms 
corresponds to the inter planar covalent bond discussed in the text. 

(dark lines in Fig. show covalent character, while the in plane electrons are delocalized. 
The density of states show a pseudo gap attributed to the covalent band. Interestingly, the 
parent structure called as Ga-II is a distorted tetragonal face centered cubic (FCC), which 
can be obtained under pressure of 35 kbar. This structure is quite close in energy to the 
cv-Ga and is metallic. In other words Ga-II, the metallic FCC structure rearranges into 
cv-Ga, once the pressure is released forming a covalent bond. Thus, the tendency to form 
a covalent bond is already there in a nascent form and given a chance to further rearrange 
the atomic positions, as is the case in clusters, the structure may prefer to form covalent 
bonds. This is what we observe in the finite sizes of Ga n clusters where, n « 55. In fact, 
our detailed analysis of the ground state geometries of Ga„ (n= 13-55) clusters shows an 
interesting pattern as to how the covalent bonds in the cluster rearrange themselves, leading 
to significant variations in the melting temperature. At what sizes the partial metallic bonds 
are formed is an interesting question but beyond the scope of this work. 

In what follows, we study the ground state geometries of selected gallium clusters (viz., 
Ga n , n= 13, 17, 20, 30, 31, 33, 37, 46, and 55) and discern the evolutionary trends. Along 
with the evolutionary trends, we also analyze the nature of bonding, the number of covalent 
bonds, and their distribution within the clusters. Using these trends, we bring out the factors 
contributing to the non-monotonic variations in the melting temperature. Thus, we expect 
our current work to throw light on general experimental observations over the reported size 
range. 
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II. COMPUTATIONAL DETAILS 



In order to have a realistic guess for the ground state, we have optimized about 300 
geometries for each of these clusters. The initial configurations for the optimization were 
obtained by carrying out a constant temperature dynamics of 100 ps each at various temper- 
atures between 600 to 1200 K. The optimization was carried out using Vanderbilt's ultra soft 
pseudo potentials^ within the Generalized Gradient Approximation (GGA) for describing 
the core-valance interactions as implemented in the vasp package.— For all the calculations, 
we take the 4s 2 and 4p l electrons as valence electrons and the 3d electrons as a part of the 
ionic core. An energy cutoff of about 9.54 Ry is used for the plane wave expansion of the 
wave function, with a convergence in the total energy of the order of 0.0001 eV. 

Once the ground state geometry is obtained from the first principles calculations, various 
structural and electronic properties of the cluster in its ground state were analyzed. The 
structural properties analyzed are (i) bond length variations within the cluster, (ii) distance 
from Centre of Mass (COM), and (iii) shape of the cluster. The shape of the cluster is 
quantified using the deformation coefficient (e pro ). For a given configuration, e pro is defined 
as 

2Q Z 



Qx Qy 

where, Q x , and Q 2 are the eigenvalues, in ascending order of the quadrapole tensor 



(1) 



Qij — RuRij (2) 
i 

/runs over the number of ions and R/j is the ith coordinate (i and j run from 1 to 3) of the 
ion 'I' relative to the center of mass of the cluster. In simple terms, Q x , Q^, and Q z define 
the spread of the cluster along the X, Y, and Z axis. Thus, configuration with spherical 
shape has eigen values Q x ~ Q y ~ Q z (e pro =l). A prolate configuration has Q 2 >> Q y ~ 
Q x , while a structure with oblate configuration has Q z ~ Q y » Q, x . 

The nature of bonding between the atoms in a cluster is analyzed using Electron Local- 
ization Function (ELF) . 1 ^ 21 For a single determinantal wave function built from Kohn-Sham 
orbitals, ipi, the ELF is defined as,™ 

Xelf = [1 + (D/D h f]-\ (3) 

where 

D h = (3/10)(37r 2 ) 5/ y/ 3 , (4) 



D = (1/2) £ |V^| 2 - (1/8)| Vp| 2 /P, (5) 

i 

with p = p(r) being the valence electron density. D is the excess local kinetic energy density 
due to Pauli repulsion and is the Thomas-Fermi kinetic energy density. The numerical 
values of Xelf are conveniently normalized to a value between zero and unity. A value 
of 1 represents a perfect localization of the valence charge while the value for the uniform 
electron gas is 0.5. Typically, the existence of an isosurface in the bonding region between 
two atoms at a high value of Xelf say, > 0.70, signifies a localized bond in that region. 

Recently, Silvi and Savin 19 introduced a nomenclature for the topological connectivity of 
the ELF. According to this description, the molecular space is partitioned into regions or 
basins of localized electron pairs. At very low values of ELF all the basins are connected 
(disynaptic basins). In other words, there is a single basin containing all the atoms. As the 
value of Xelf is increased, the basins begin to split, and finally we will have as many basins 
as the number of atoms. The value of ELF at which the basins split (a disynaptic basin 
splits into two monosynaptic basins) is a measure of the interaction between the different 
basins (a measure of the electron derealization). 

III. RESULTS AND DISCUSSION 

In order to gain some insight into the experimentally observed differences in the heat 
capacity curves of Ga n (n = 17 — 55), we have investigated the ground state geometries 
of selected clusters, viz., Ga n , n = 13,17,20,30,31,33,37,40,44,46,55 and analyzed the 
differences in their structure and bonding. The choice of these clusters is dictated so as to 
represent the changing nature of the heat capacities across the measured series. For example, 
among these clusters, Ga^, Ga3o, and Gass have very broad heat capacities. Ga4o and Ga44 
show a weak peak around 550 K and 700 K, respectively. Ga2o, Ga3i, Ga33, Ga37, and Ga46 
have a well recognized peak. Among these, Ga2o and Ga46 melt around 800 K, while Ga3i, 
Ga33, and Ga37 melt between 550-600 K. Gai3 is taken as a reference cluster to analyze the 
growth pattern in these clusters. 

In Fig. |2l we show the ground state geometry of these clusters obtained in our search. The 
atoms shown in light color (blue on line) indicate the growth of the cluster over the previous 
one. The shape of these clusters is analyzed using e pro and eigen values of quadrapole tensor 
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(a) Ga i3 (b) Ga u 



(c) Ga 20 (d) Ga 30 



(e) Ga 31 




(f) Ga 3 i 



(g) Ga 33 



(h) Ga 37 



(i) Ga 40 




(j) Ga 



■ ii 



(k) Ga 46 



(1) Ga, 



55 



FIG. 2: Ground state geometries of selected gallium clusters between Gai 3 and 



which are shown in Fig. |3] and Fig. respectively. We begin our discussion with a note 
on Gai 3 shown in Fig. [2]-(a) . The ground state geometry of Gai 3 was found to be slightly 
distorted decahedron in previous works . 10 i 13 i 20 The cluster is nearly spherical in shape with 
1.1. The inter-planar and core-surface bonds are the shortest ones (2.43-2.65 A) 111 



'-pro 



the cluster. The intra-planar bond distances are the longer ones and vary between 2.85- 
3.00 A. Addition of 4 atoms to Gai 3 results in a deformation of decahedron as seen from 
Fig. EJ-(b). The atoms shown in red (dark color) indicate the deformed structure of Gai 3 , 
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FIG. 3: Deformation co-efficient parameter (e pro ) of various gallium clusters. 



while the atoms in blue (light color) show the extra four atoms. It is clearly seen from this 
figure as well as from the analysis of Q x , Q y , and Q 2 (see Fig. Hj) that Gai7 is a growth over 
Gai3 along Z-axis. The e pro of Gai7 is ~ 2.1. The shortest bonds (2.50-2.70 A) in Gai7 are 
randomly distributed (scattered) within the cluster unlike in case of Gai3, where they are 
localized. The implication of this distribution of shortest bonds is brought out more clearly 
in the later section where we discuss their connectivity and its relation with the melting 
temperature. 

Fig. EHc) shows the ground state geometry of Ga2o- The atoms in light color (blue on 
line) are the extra seven atoms that form a dome like structure on the deformed Gai3 cluster. 
The cap to the top plane of Gai3 now becomes the core atom in case of the Ga2o- Thus, 
Ga2o is also obtained by addition of atoms on deformed Gai3 along the Z-axis. The eigen 
values of quadrapole tensor and e pro (~ 1.8) clearly indicate Ga 20 to be prolate in shape. 
As in case of Gai3 and Gai7, the atoms are bonded to each other through bond lengths 
ranging between 2.54-3.00 A. A careful examination of the distribution of bond lengths in 
Ga2o cluster reveals that the shortest bond lengths (2.54-2.62 A) in the cluster are localized 
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Size of the Cluster 



FIG. 4: Eigen values of the quadrapole tensor for gallium clusters. 

in the newly added (top) and bottom planes of the cluster. The core atom of the cluster is 
connected to the surface atoms with bond lengths of 2.65-2.70 A, while rest of the atoms 
within the cluster are connected to each other through varying lengths of 2.73-3.00 A. 

Addition of 10 atoms to Ga2o, viz., Ga3o, results in a configuration with stacked planes as 
shown in Fig.El-(d). The Q x , Q y , and Q z values show that the growth is predominantly along 
the Y and Z-axis resulting in a nearly oblate configuration. Ga3o does not have a distinct 
core. However, it has two interstitial atoms (atoms not facing the surface) connected to each 
other and to the surface atoms through bond lengths of ~ 3.0 A (shown in light red color on 
line). The shortest bonds are randomly distributed within the cluster. Addition of one atom 
to Ga3o results in a significant reordering of the planes and symmetric distribution of atoms 
in the cluster (cluster has Co v symmetry) as seen from Fig. 0-(e) (for a detailed discussion 
on Ga 30 and Ga 3i see Ref.llll). Another perspective of Ga 3 i (Fig. I2-(f ) ) and Q x , Q y , and 
values clearly show the cluster to be oblate in shape. Ga 3 i also does not have a distinct core. 
It has three interstitial atoms (shown in light red color on line) which are weakly bonded 
(distances > 2.90 A) to each other as well as to surface atoms. The shortest bonds in the 
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FIG. 5: Distance of atoms from Centre of Mass (COM) in Ga3i and Ga46 clusters. 

cluster are localized within the same plane. Fig. Eh(g) shows Ga33 along the X-Z plane. 
We note that while Ga33 remains identical to Ga3i along the Y-Z plane (figure not shown), 
what is visibly different is the presence of an extra protruding atom along the X-axis. This 
growth continues in Ga37 as seen from its ground state geometry in Fig. EJ-(h). This growth 
is also evident from the increasing values of Q x and a sudden decrease in the e pro . Both 
Ga33 and Ga37 retain the stacked planes configuration similar to Ga3o and Ga3i. The planes 
are well defined in both the cases and the clusters do not have a distinct core. The clusters 
continue to grow along the the X-axis as seen in the ground state geometry of Ga4o (shown 
along the X-Z plane) in Fig. I2h(i)- Addition of four more atoms to Ga4o, however, leads 
to a significant rearrangement of atoms within the cluster with out any preferential growth 
along any axis resulting in a nearly spherical structure of Ga44 (shown in Fig. EHj)) with e p 
~ 1.1. We note that as the cluster evolves into a spherical shape as in case of Ga44, we once 
again note the presence of a distinct core shell. Ga44 has five core atoms (shown in light 
red color on line). Addition of 2 atoms to Ga44 results in a highly symmetrical cluster (C3„ 
symmetry) with 9 core atoms and a spherical configuration (Ga46 shown in Fig. EHk)) with 
e p ~ 1.0. To summarize, the clusters undergo a spherical to prolate transition (Gai3-Ga2o) 
followed by a prolate to oblate transition (Ga 20 -Ga 3 i) and an oblate to spherical transition 
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(Ga3i-Ga46). As the clusters grow further, it is seen that the atoms begin to accumulate 
once again preferentially along one axis, leading to prolate configurations as seen in Gas5 
(Fig. [2H1))- We predict here that the clusters are likely to undergo once again a prolate to 
oblate and an oblate to spherical transition as the number of atoms continue to increase. In 
the studied size range we notice that the clusters with spherical configurations have distinct 
core, while the clusters with oblate configurations do not have any core. This is clearly seen 
from Fig. in which we show the distance of atoms from the COM of Ga3i which is oblate 
in shape and Ga 46 , a spherical structure. Ga 46 has a clear shell structure, while Ga 3i does 
not have a distinct core but is characterized by interstitial atoms. 

To understand the reasons behind the occurrence of so called "melters" and "non- 
melters", we analyze the nature of bonding in the clusters studied. In Fig. El we show 
the isosurface of ELF at a value of Xelf = 0.70 for all the clusters. The black lines join the 
atoms in the same basin. A careful examination of the basin structures reveals an interesting 
contrast between "non-melters" (Ga n , n= 17, 30, 40, 44, and 55) and "melters" (Ga n , n=20, 
31, and 37). The "non-melters" show a fragmented pattern of basins. In other words, there 
are many basins each consisting of not more that 4-5 atoms. In contrast, "melters" have at 
least one large basin consisting of 10 atoms or more. Ga46, a "melter", apparently appears 
to have fragmented basin as seen visually. A detailed examination reveals that all the 9 core 
atoms are connected at the value of Xelf = 0.70 forming a single basin (indicated by an 
arrow). We further note that these core atoms are also connected to the surface atoms, there 
by forming a single large basin of 15 atoms. This is made more clear in the next section 
where we show the schematic representation of the covalent bonds in each cluster. 

To summarize, the clusters exhibiting a sharp solid-liquid transition are seen to have 
substantial number of atoms connected through a single basin. Thus, a large group of 
atoms are bonded together with a similar strength forming an island of local order and 
therefore it is reasonable to expect that they will 'melt' together. In this sense the cluster 
can be considered as (at least partially) ordered and will show a well defined peak in the heat 
capacity. On the contrary, in clusters exhibiting broad or continuous solid-liquid transition, 
each atom (possibly a group of atoms) has different local environment. That means different 
atoms are bonded with the rest of the system with varying strength. Consequently, their 
dynamical behavior as a response to temperature will differ. Some of the atoms pickup 
kinetic energy at low temperatures, while the others may do so at higher temperatures. 
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Gaii <T m ~ 600) Gaai (T m ~ 550) Ga^ (T m ~ 600) 




Gas (T m ~ 800) Gaje (T m ~ 800) 

FIG. 7: Distribution of shortest bonds within the gallium clusters. 

Therefore, we can expect the cluster to have a broad and continuous melting transition. 
Thus, our results confirm our earlier argument that there exists a definitive relationship 
between the local order in the cluster and its finite temperature behavior (and consequently 
the characteristics of the heat-capacity curve).— 

We also show that as the cluster grows in size, it evolves through a succession of ordered 
and disordered geometries. In such cases an addition of one or few atoms changes the 
nature of the ground state geometry abruptly. The variation in the ground state results in 
the appearance or disappearance of the local order within the structure, leading to presence 
or absence of the melting peak, respectively. We now realize that this size sensitive nature 
of heat capacities is generic to small clusters and related to the evolutionary pattern seen 
in their ground states. The evidence for this comes not only from gallium clusters but also 
from clusters of sodium , 14 ' 2 ? aluminum,— and gold^ having very different nature of bonding. 
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FIG. 8: Variation in the binding energy (de) of cluster with respect to the displacement of surface 
atoms from their equilibrium positions (Ro). 

Next, we attempt to understand the reason behind the variations in the melting temper- 
ature in Ga n (n = 20-46) clusters. We analyze the structure and bonding of the clusters 
exhibiting a distinct melting transition viz., Gai3, Ga2o, G&31, Ga33, Ga37, and Ga46. It 
may be recalled that Ga2o and Ga46 melt around 800 K, while Gagi, Ga33, and Ga37 melt 
between 500-600 K. The first striking difference to note is the presence of core atoms in the 
Ga2o and Ga46 and absence of any distinct core in Ga3i, Ga33, and Ga37. The distribution 
of the shortest bonds (< 2.70 A) in each cluster (which are also the atoms forming a single 
basin as indicated by ELF) is shown in Fig. [7| It is clearly noted that the shortest bonds 
in Ga2o and Ga46 are spread through out the cluster and more importantly they make a 
network connecting the core and the surface atoms. This results in a compact and stable 
electronic configuration in both the clusters as indicated by bunched eigen value spectrum 
as compared to evenly spread spectrum in case of all other clusters (figure not shown). On 
the other hand, in case of Ga3i, Ga33, and Ga37 the covalent bonds are restricted along 
the planes of the cluster. Thus, the pattern of bonding within these three clusters appears 
similar to that of pillared materials, where the atoms along the same plane are covalently 
bonded, while the bonds across the planes are weakly covalent or more metallic like in na- 
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ture. Hence, when these clusters are heated, the weaker inter-planar bonds are the first ones 
to break. This is followed by sliding of planes (constituting of equivalently bonded atoms) 
along each other. This argument correlates well with our analysis of ionic motion in Ga3i,— 
where we see that as the cluster is heated, the planes begin to slide as shown by arrows 
in Fig. [7| This is initiated at a fairly lower temperature of 400 K. On further increase in 
the temperature of the cluster, the sliding planes collide against each other and the cluster 
melts around 600 K. On the other hand, the strong bonds between the core and surface 
atoms in Ga2o hold the surface atoms around their respective positions until a temperature 
of 600 K. 12 The cluster melts finally around 800 K. We further note that another cluster 
with a distinct core viz., Gai3 also exhibits a network of strong covalent bonds between the 
core and surface atoms. Our earlier simulations indicate this cluster to melt at a relatively 
high temperature of 1300 K.— 

Thus, we believe that the presence of a network of strong covalent bonds between the 
core and surface atoms results in greater stabilization of the surface atoms in Ga2o and 
Ga46- These surface atoms therefore need to overcome a higher energy barrier to diffuse 
across the surface as compared to one in Ga3i, Ga33, and Ga37. To demonstrate this, we 
stretch one of the atoms on the surface of Ga3i and Ga46. This is done so as to mimic the 
motion of the atoms when the cluster is heated. The change in the binding energy of the 
cluster as function of the bond stretch is plotted in Fig. [7| It is easily noted from the figure 
that the surface atom of Ga46 requires more energy for displacement from its equilibrium 
position as compared to the one in Ga3i. This is true for most of the surface atoms of 
Ga46 and Ga3i. This analysis is also found to be similar in case of aluminum and sodium 
clusters. Our understanding is also in agreement with the results of Aguado and Lopea^ for 
sodium clusters, where they have found a strong correlation between variation in T m and 
core-surface distances. 

IV. SUMMARY AND CONCLUSION 

In the present work we have attempted to explain the reasons behind the characteristic 
features observed in the experimental heat capacity curves of gallium clusters in the size 
range of 20-55. As the gallium clusters grow, they are seen to pass through a cycle of 
spherical-prolate, prolate-oblate, and oblate-spherical transition. During the process they 
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pass through a succession of geometries with and without local "order" in bonding. It is the 
presence of a local "order" in the cluster that is responsible for the sharp peak in the heat 
capacity curve. Our studies also show that it is the presence of a distinct core and more 
importantly presence of a network of strong bonds between the core and surface atoms that 
is responsible for the higher melting temperature of the clusters. 
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